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q We report the discovery of the two-image gravitationally lensed quasar SDSS J1335+0118. The object 

was selected as a lens candidate from the Sloan Digital Sky Survey. The imaging and spectroscopic follow- 
ed , up observations confirm that the system exhibits two gravitationally lensed images of a quasar at z = 1.57. 

The image separation is 1"56. We also detect an extended component between the two quasar images, 
likely the lensing galaxy. Preliminary mass modeling predicts the differential time delay At ~ 30/i _1 day 
assuming the redshift of the lens galaxy is 0.5. 
/\ ' Key words: gravitational lensing — quasars: individual (SDSS J133534. 79+01 1805.5) 

1. Introduction quasars can be a powerful tool to study quasar it- 
self through microlensing, spectra, or variabilities of 

Multiple images of a quasar produced by a foreground lensed quasars (Chang & Refsdal 1984; Grieger, Kayser, 

galaxy are useful for cosmological and astrophysical ap- & Schramm 1991; Mediavilla et al. 1998; Mineshige 

plications, including the measurement of the cosmological & Yonehara 1999; Yonehara 1999; Wyithe, Webster, 

constant through the lensing rate (Turner 1990; Fukugita, & Turner 2000; Yonehara 2001; Lewis & Ibata 2003). 

Futamase, & Kasai 1990; Chae et al. 2002), the determina- Therefore, it is important to make a large sample of lensed 

tion of global Hubble constant from time delays between quasars also to find more gravitational lens systems which 

images (Refsdal 1964, 1966), the study of the mass distri- are suitable for studying quasar. 

butions of the lensing galaxies (Kochanek 1991; Metcalf The Sloan Digital Sky Survey (SDSS; York et al. 2000; 

& Madau 2001; Chiba 2002; Dalai & Kochanek 2002; Stoughton et al. 2002; Abazajian et al. 2003) has the 

Schechter & Wambsganss 2002; Koopmans & Treu 2003; potential to provide such a lensed quasar sample. The 

Rusin, Kochanek, & Keeton 2003), and the formation and SDSS is a survey to image 10 4 deg 2 of the sky as well as 

evolution of galaxies (Kochanek et al. 2000; Kochanek to obtain spectra of galaxies and quasars from the imag- 

& White 2001; Keeton 2001a; Oguri 2002; Ofek, Rix, ing data. The dedicated 2.5 meter telescope at Apache 

& Maoz 2003). A large homogeneous sample of lensed Point Observatory is equipped with a multi-CCD camera 

quasars with a well-understood selection function is cs- (Gunn et al. 1998) with five broad band centered at 3561, 

sential in these statistical studies. In addition, lensed 4676, 6176, 7494, and 8873A (Fukugita et al. 1996). The 
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imaging data are automatically reduced by the photomet- 
ric pipeline (Lupton et al. 2001). The astrometric posi- 
tions are accurate to about O'.'l for sources brighter than 
r = 20.5 (Pier et al. 2003). The photometric errors are typ- 
ically less than 0.03 magnitude (Hogg et al. 2001; Smith 
et al. 2002). The SDSS quasar selection algorithm is 
presented in Richards et al. (2002); the SDSS spectra 
cover 3800-9200 A at a resolution of 1800-2100. The fi- 
nal spectroscopic quasar sample is expected to comprise 
10 5 quasars, thus will contain ~ 10 2 lensed quasars given 
the typical lensing probability 0.1% (Turner, Ostriker, & 
Gott 1984). Indeed, several new gravitational lens sys- 
tems have been discovered using SDSS data (Inada et al. 
2003a,b,2004; Morgan, Snyder, & Reens 2003; Johnston 
et al. 2003; Pindor et al. 2004; Oguri et al. 2004). 

We report the discovery of the lensed quasar 
SDSS J133534.79+011805.5 (SDSS J1335+0118) in the 
SDSS. This quasar has already been identified in the 
Large Bright Quasar Survey (LBQS 1333+0133; Hewett 
et al. 1991), but not as lensed. A search for lens candidates 
in the LBQS was carried out by Hewett et al. (1998), but 
those authors state that their search is not sensitive to lens 
systems with < 3" separations. The quasar was selected 
as a lens candidate in the course of an ongoing search 
for strongly lensed quasars in the SDSS. From the results 
of photometric and spectroscopic follow-up observations, 
we conclude that the quasar is lensed by an intervening 
galaxy. 

2. Observations 

2.1. Candidate Selection 

To select SDSS lensed quasar candidates, we examine all 
spcctroscopically confirmed quasars with z > 0.6. Quasars 
at smaller redshifts have much lower probability of being 
lensed and also are often extended because of their host 
galaxies, making it more difficult to select lens candidates 
with our candidate selection algorithm described below. 
We search for lens candidates using a combination of SDSS 
parameters: dev_L (the likelihood parameter of fitting by 
dc Vaucoulcurs profile), exp_L (the likelihood parameter 
of fitting by exponential disk), and star_L (the likelihood 
parameter of fitting by point spread function). This can- 
didate selection algorithm has already found three new 
lensed quasar systems (Inada et al. 2003a, 2004; Pindor et 
al. 2004), and can identify gravitationally lensed quasars 
with image separations lf'0 — 2"5 quite well (N. Inada et 
al., in preparation). The SDSS ugriz imaging data of 
SDSS J1335+0118 (sec Figure 1) clearly shows that the 
system consists of two stellar components with similar col- 
ors, making it an excellent lensing candidate. The total 
magnitudes of this system in the SDSS photometric data 
are u = 18.21 ± 0.02, g = 17.83 ± 0.01, r = 17.62 ± 0.01, 
i = 17.26 ±0.01, and z = 17.14 ±0.03. 

2.2. Additional Imaging Observations 

We obtained a deep i-band image on 2003 May 28 with 
the Subaru Prime Focus Camera (Suprime-Cam; Miyazaki 
et al. 2002) of the National Astronomical Observatory of 




Fig. 2. Subaru Suprime-cam images of SDSS J1335+0118 in 
the i-band. The image scale is 0''2pixel — 1 . Left: original 
Suprime-cam image. Components A and B are lensed quasar 
components. Right: PSF-subtracted image. The component 
G is likely to be the lens galaxy. 

Japan's Subaru 8.2-m telescope, which is shown in Figure 
2. The exposure time was 30 seconds and the seeing was 
0'.'5 — 0'.'6. The image scale is 0'.'2pixel _1 . Each frame 
was bias-subtracted and flat-field corrected. Subtraction 
of the two quasar images using a nearby star as a point- 
spread function (PSF) template reveals an extended ob- 
ject, denoted as G, located on the line between A and B 
and closer to B (see Figure 2). This configuration is ex- 
pected for a standard simple lens model such as singular 
isothermal sphere model, so component G is likely to be 
the lens galaxy. In Table 1, we summarize the position 
of each component and the relative fluxes of the quasar 
components that are measured using Subaru Suprime-cam 
data. The errors of the galaxy position is determined from 
Gaussian fits. 

We also obtained J and -fT-band data on 2003 April 
19 with the Near InfraRed Camera (NIRC; Matthews & 
Soifer 1994) of the Keck I telescope at the W. M. Keck 
Observatory on Mauna Kea, Hawaii, USA. The total ex- 
posure time was 900 seconds, seeing was variable, av- 
eraging ~ 0'.'6. Conditions were not photometric, with 
about 1 magnitude of variable cloud. The image scale 
is 0''15pixel _1 . Standard reduction procedures were fol- 
lowed to remove the dark current and to flat field the data. 
We show the K image in Figure 3. Although the galaxy 
component can be seen in the direct image, we subtracted 
the PSF of two quasar components by fitting a Gaussian 
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Table 1. Positions and flux ratios of SDSS J1335+0118 in the Subaru i-band image 



Object 



y [arcsec] * Flux [arbitrary] 



A 0.000 ±0.001 0.000 ±0.001 1.0 ±0.2 

B -1.038 ±0.002 -1.165 ±0.002 0.374 ±0.075 
G -0.769 ±0.011 -0.757±0.011 

* The positive directions oi x and y are defined by West and North, respectively. 
t Errors are broadened to 20% to account for possible systematic effects. 



K(subtracted)" 




Fig. 3. Same as Figure 2, but the J<"-band images taken with 
NIRC on Keck I telescope are shown. The image scale is 
0."15pixel -1 . 

PSF at each quasar location, revealing the lensing galaxy 
more clearly (see Figure 3). As in the i-band image shown 
in Figure 2, the galaxy component is seen in the PSF- 
subtracted if -band image. We compare the distance be- 
tween components A, B, and G in Table 2. The differences 
of the distances are large, but only ~0."1 at most and still 
much smaller than seeing sizes and pixel sizes. The po- 
sition of the lens galaxy is particularly uncertain in the 
K data because not only is the lensing galaxy relatively 
faint and extended, but also because the seeing was quite 
variable on the night the IR data were obtained and no 
suitable PSF star is present in the NIRC field, making it 
difficult to properly subtract the PSF wings which over- 
lap the galaxy. Thus it may be possible that such residu- 
als systematically affect the estimated positions. We note 
that the system is detected (but not resolved) by the Two- 
Micron All-Sky Survey (2MASS) with J = 16.23 ±0.10, 
H = 15.47 ±0.11, and K = 15.29 ± 0.18. 

2.3. Spectroscopic Observations 

The spectroscopic observation was done on 2003 March 
2 with the red CCD of the ESO Multi-Mode Instrument 
(EMMI) on the ESO New Technology Telescope (NTT). 
We used grism #3 (3601inemm _1 , 2.3Apixcl~ 1 , cover- 
ing 3850 - 8450 A). The seeing was l'/O. The exposure 
time was 900 sec. The slit was positioned to obtain spec- 
tra of the two quasar components (A and B) simultane- 
ously. The two components were clearly separated in the 
EMMI image, making extraction simple. The data were 
reduced in a standard way using IRAF 1 . The spectra are 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



shown in Figure 4. We find that both components A 
and B have ClV, Cm], and Mgll emission lines at the 
same wavelengths. The redshifts of components A and 
B estimated from Mgll emission lines are z = 1.57 ±0.03 
and z = 1.57 ±0.05, respectively. We cross-correlate the 
two spectra and estimate the velocity difference between 
two quasar components to be 20±800kms~ 1 . We also 
plot the ratio of the fluxes of components A and B, and 
find that the ratio is almost constant (~ 0.35) for a wide 
range of wavelengths. Therefore we conclude that com- 
ponents A and B are two gravitationally lensed images of 
a quasar at z = 1.57. The lens hypothesis is further sup- 
ported by the fact that both components appear to have 
associated ClV absorption at a blueshift of 2500kms _1 
from the quasar. We note that both A and B have strong 
Fell/Mgll/Mgl absorption system at z = 1.43, Alii ab- 
sorption system at z = 1.42. This absorption is unlikely 
to be associated with the lens, given the small difference 
between the absorption and source redshifts. Indeed, the 
lensing galaxy is typically half way out in affinc or an- 
gular diameter distance, and probabilities of such small 
difference between source and lens redshifts are too small 
(e.g., Ofek, Rix, & Maoz 2003). Moreover, such absorp- 
tion is quite common since the typical linear separation 
between quasar sightlines and galaxies which are respon- 
sible for such absorption is large, ~ 50kpc (Bergeron & 
Boisse 1991; Fukugita, Hogan, & Peebles 1996). It seems 
that strengths of absorptions are different between com- 
ponents A and B, but this is not surprising given the in- 
homogeneity of the spatial distribution of absorbers (e.g., 
Petitjean & Bergeron 1990; Steidel & Sargent 1992). We 
note that the color of component G (see §2.4) is consis- 
tent with that of a low-redshift (z < 0.5) early-type galaxy 
(e.g., McLeod et al. 1995). 

2-4- Differential Extinction 

We check colors of quasar components to see whether 
the system suffer from differential extinction. The results 
of photometry from the imaging data of the SDSS, Subaru, 
and Keck are shown in Table 3. Relative magnitudes are 
estimated from Subaru, Keck, and SDSS images. To cal- 
ibrate total magnitudes, we use the measurements of the 
SDSS, i = 17.26 ±0.01, and the 2MASS, K = 15.29 ±0.18. 
We find that the colors of each component are consistent 
with a single value, given the errors of ~ O.lmag associ- 
ated with the deconvolution of two components. Spectra 
of quasar components shown in Figure 4 also indicates 
no significant differential extinction. Although it seems 
that at < 4800 A the ratio increases as the wavelength de- 
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Table 2. Comparison of positions between Subaru i-band and Keck K -band images 



Objects i-band Distance [arcsec] Jv-band Distance [arcsec] 
AB 1.560 ±0.002 1.503 ±0.002 

AG 1.079 ±0.011 0.969 ±0.012 

BG 0.489 ±0.011 0.542 ±0.012 



2000 3000 




4000 5000 6000 



Fig. 4. ESO NTT spectra of SDSS J1335+0118 components 
A and B. The redshifts of components A and B estimated from 
Mgn emission lines are z = 1.57±0.03 and z = 1.57±0.05, re- 
spectively. Both components have associated C IV absorption 
at a blueshift of 2500km s — 1 from the quasar (marked by ©). 
In addition, both components have strong Fe Il/Mg Il/Mg I 
absorption system at z = 1.43, and Alii absorption system at 
z = 1.42. Note that the absorption at ~ 7600A is atmospheric. 
The bottom panel shows the ratio of spectra. 

creases, but the feature is not significant because of the 
large errors in those wavelengths arising from the low re- 
sponse at ~ 4000 A or shorter. Therefore we conclude that 
differential reddening is not significant. 

3. Mass Modeling 

We model the lens system with the simple Singular 
Isothermal Ellipsoid (SIE) model, assuming that compo- 
nent G is the lens galaxy. The SIE model has the following 
dimcnsionlcss surface mass density: 

*(r) = ^, (1) 

where ( = r[l + ((l-q 2 )/(l + <f))cos2(e-6 e )} 1 / 2 , r cin is 
the Einstein ring radius, q is the lens axis ratio, and 9 e is 
the position angle of the ellipse. The model includes eight 
parameters: the galaxy position x g and y s , the Einstein 
ring radius r e j n , the ellipticity e = 1 — q, the position an- 
gle 6 e , the source position x s and y a , and the flux of the 
quasar f 2 . Since the number of constraints is also eight 

2 We obtain the same results even if we do not adopt a parameter 
/ and instead we use only the flux ratio, not fluxes of A and B, 
as a constraint. 



Table 4. Best fit SIE model parameters for 
SDSS J1335+0118 



Parameter 


Value 


r ; n [arcsec] 
e 

# e [deg] 


u - ' oy -0.011 

n 1 q+o.04i 

U.1ZO_q 024 
n n + 19.3 
_ '-°-12.0 



Wc use positions and flux ra- 



tios of SDSS J1335+0118 in 
the Subaru i-band image (see 
Table 1). Errors are 68% con- 
fidence. 

(see Table 1), the degree of freedom is zero. We use the 
public software lensmodel (Kecton 2001b) to constrain 
the models. The result is summarized in Table 4. The 
best-fit parameters predict the time delay between im- 
ages as At = 29.5ft.~ 1 day (h denotes the Hubble constant 
in units of 100km s~ 1 Mpc -1 ) if we adopt a lens redshift 
0.5, flM = 0.3, and Oa = 0.7. Errors on the parameters are 
estimated from one-dimensional slices of the y 2 surface. 
We find that we can fit the lens system with relatively 
small lens ellipticity e = 0.123l£^, significantly smaller 
than observed ellipticity of the light (~ 0.3). The derived 
position angle of the lens galaxy is consistent with that 
of the light (~ 7°) within la. These results are in good 
agreement with previous studies (e.g., Keeton, Kochanek, 
& Falco 1998). 

4. Conclusion 

We have reported the discovery of a new two-image 
gravitational lens SDSS J1335+0118. The system was 
identified as a new lens candidate in the SDSS. The im- 
age separation is 1"56± 0.002 in the i-band image taken 
with Subaru Suprime-cam. The spectroscopic observation 
with ESO NTT has confirmed both images have identical 
redshift z = 1.57 (z = 1.57 ±0.03 and z = 1.57 ±0.05 for 
components A and B, respectively). We have also proba- 
bly identified the lensing galaxy between two quasar com- 
ponents in the higher-resolution imaging data taken with 
Subaru Suprime-cam and Keck NIRC. The lens geometry 
is well reproduced with a simple mass model and reason- 
able model parameters. Assuming the redshift of the lens 
galaxy as 0.5, we have predicted a differential time delay 
At ~ 30/i- 1 day. 

We thank an anonymous referee for many useful com- 
ments. Funding for the creation and distribution of the 
SDSS Archive has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National 
Aeronautics and Space Administration, the National 
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Table 3. Photometry of Subaru i-band, Keck _R"-band, and SDSS ugriz-band images 



Band 


Total 


A 


B 


B-A 


G 


u (SDSS) 


18.21 ±0.02* 


18.54 


19.68 


1.14 




9 (SDSS) 


17.83 ±0.01* 


18.12 


19.39 


1.27 




r (SDSS) 


17.62±0.01* 


17.95 


19.09 


1.14 




i (SDSS) 


17.26±0.01* 


17.60 


18.68 


1.08 




i (Subaru) 


17.26±0.01* 


17.63 


18.70 


1.07 


20.05 


z (SDSS) 


17.14±0.03* 


17.49 


18.54 


1.05 




K (Keck) 


15.29±0.18 t 


15.78 


16.75 


0.97 


17.81 


Based on the SDSS. 
Based on the 2MASS. 
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